The Equations of Motion of a Charged Particle in the 
Five-Dimensional Model of the General Relativity Theory with 
the Four-Dimensional Nonholonomic Velocity Space 



We consider the four-dimensional nonholonomic distribution defined by the 4- 
potential of the electromagnetic field on the manifold. This distribution has a metric 
tensor with the Lorentzian signature (+, — , — , — ), therefore, the causal structure ap- 
pears as in the general relativity theory. By means of the Pontryagin's maximum 
principle we proved that the equations of the horizontal geodesies for this distribution 
are the same as the equations of motion of a charged particle in the general relativity 
theory. This is a Kaluza - Klein problem of classical and quantum physics solved by 
methods of sub-Lorentzian geometry. We study the geodesies sphere which appears in 
a constant magnetic field and its singular points. Sufficiently long geodesies are not 
optimal solutions of the variational problem and define the nonholonomic wavefront. 
This wavefront is limited by a convex elliptic cone. We also study variational princi- 
ple approach to the problem. The Euler - Lagrange equations are the same as those 
obtained by the Pontryagin's maximum principle if the restriction of the metric tensor 
on the distribution is the same. 
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1 Nonholonomic model in sub-Lorentzian geometry 

Distribution on a smooth manifold M is a family of subspaces A(x) C T X M (of the tangent 
bundle of the manifold) with the same dimension, smoothly parametrized by the points of 
manifold [U [21 El H]. For each x G M on A(x) the bilinear form (u,u) x is defined which 
can be equally described by the metric tensor gij(x). The metric tensor of a distribution 
smoothly depends of the point of the manifold. The absolutely continous path x(t) is called 
horizontal, iff x'(t) G A(x(t)) for almost any t. In sub-Riemannian geometry the metric 
tensor gij(x) of the distribution A is positively defined. In this geometry we study the 
shortest horizontal curves connecting two given points. If the distribution is completely 
nonholonomic and the Riemannian manifold is full, then any two points of the manifold can 
be connected by the shortest horizontal geodesic. On manifolds and also on distributions 
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with metric tensor the geodesic is the Euler - Lagrange solution for the length functional 

J(x{-)M-)) = !lMt)Mt))%dt. 

In sub-Lorentzian geometry the metric tensor gij{x) of the distribution A has the signa- 
ture (+, — )• In this geometry we study the longest horizontal curves connecting two 
given points |5J. The equations of motion of a charged particle in the electromagnetic and 
gravitational fields can be obtained as the solutions of the variational problem for the dis- 
tribution |6j. Consider the 5-dimensional manifold M 5 with the 4-dimensional distribution 

3 

A defined by means of the differential form uj(x) = Ai(x)dx l + dx 4 , where the covector 

i=0 

(ylj)j =0 3 is the 4-potential of the electromagnetic field. Consider the length functional 

J(x(-), «(•))= / L{x{t),u{t))dt, (1) 



J t 

1 /2 

where L(x, u) = m(u, u) x is the pseudonorm of the vector u, m is the mass of the particle. 
Consider the problem of maximizing the length functional on the set of horizontal curves 
connecting two given points. Solution of this problem as we proove in this paper leads to 
equations 

d dL dL v-^ , 

iS-S+f'EV^, * = 0,..,3, (2) 

j=0 

where Fjk are defined by (|5|. These equations are the same as the equations of motion of 
a particle with the charge in the electromagnetic and gravitational fields of the general 
relativity theory. This wording of the problem is different from the classical problem of 
physics [7J [8] of minimizing the action functional on a 4-dimensional manifold: 

S(x(-),u(-)) = -mc rMt),u(t)>igdfc-- ry^A t {x{t))u\t)dt, (3) 
Jo c J , =o 

where e is the charge of the particle, c is the speed of light. This paper is devoted to the 
problem of unification of electromagnetic, gravitational and other fundamental interactions 

In this model the particles are considered on the 5-dimensional manifold M 5 with the 
special smooth structure. Allowed are only smooth coordinate transformations for which 

*_0, * = <),.. .,3, g-l. (4) 

Components q^j, i = 0, . . .,3, can be interpreted as gauge transformations (j9|. Since the 

coordinate transformations are smooth, Tp-j^r^ = i^-r^r-j = 0, i, j = 0, ...,4. There- 
ax 4 ox 3 ox 3 ax 4 

fore the cylindricity condition = for any tensor field v is invariant. In the tangent 

bundle TM the subspace Lin {c^} is invariant. In the cotangent bundle T*M the subspace 
Lin {dx°, dx 1 , dx 2 , dx 3 } is invariant. Assume that the metric tensor and the 4-potential of 
the electromagnetic field do not depend of x 4 . This condition is invariant at the coordinate 
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transformations with the property Q. The 4-dimensional distribution A is the velocity 
space, i.e. the set of all possible velocity vectors of the particle. The metric tensor of the 
distribution A defines 4-dimensional convex elliptic cone, completely similar to the cone 
of future of the general relativity theory. Causal structure in the proposed theory is simi- 
lar to the causal structure of the general relativity theory (T2J [13] . In the presence of the 
electromagnetic field the manifold A is nonholonomic. 

The electromagnetic field is defined in physics by the antisymmetric tensor (-Fjfc)j,fc=o,..,3 



ft 



jk)j,k=Q,..,3 



/ 
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(5) 



Here E is the 3-dimensional vector of the electric field tension, H is the 3-dimensional vector 
of the magnetic field tension. There is the 4-potential (vli)i=o,..,3 which defines the tensor of 
the electromagnetic field: 

dA k . 8Aa , . 

4-Potential is not defined unambiguously. It is defined with an arbitrary gauge transforma- 

Qf 

tion A{ t— > Ai + Qji where / is any smooth function. 

Let U be some coordinate neighbourhood on the smooth manifold M 5 . The 4-dimensional 
distribution A on U can be defined by the differential form 

3 

Ux = ^2 A i (x)dx i + dx A . (7) 

i=0 

This definition is correct since at the coordinate transformation with the property Q we 
have 

Adx* + dx 4 = ( E + ttH w + d y 4 = tl x ^y J + d y 4 - ( 8 ) 

i=0 j=0 \i=0 y y J j=0 

Therefore the coordinate transformations on the 5-manifold lead to transformations of the 
4-potential 

that include 4-dimensional coordinate transformations of the general relativity theory and 
gauge transformations of the 4-potential of the electromagnetic field [0] as in [7]. Locally 

the distribution A can be defined by the basis vector fields = -^-r — Attt, % = 0, . . ., 3. 

ox ox 

Vector fields on a manifold can be considered as operators of differentiation. During the 
coordinate transformation with the property Q we have 
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where ej = -^-j — Aj-^^, j = 0, . . ., 3. Therefore the 4-dimensional basis of the distribution 

A is transformed (10) exactly the same way as the coordinate vector fields m 
the general relativity theory. The commutators of vector fields e^, % = 0, . . ., 3, generate the 
tensor of the electromagnetic field: [e«, ej] = ~^i~^A- ^ ^ 0> then the distribution A is 
completely nonholonomic. Allowed curves on the distribution are absolutely continuous and 
satisfy almost everywhere the horizontality condition 

u 7( t)(Y(*)) = o. (ii) 

In this paper we consider the optimization problem for the length functional on the 5- 
dimensional manifold with the 4-dimensional distribution A. We prove that the solutions of 
the variational problem for this distribution satisfy the equations of motion for the charged 
particle of the general relativity theory. We consider also the geodesic sphere for the distri- 
bution A. The geodesic sphere with the radius r and center x is the set of points at distance 
r from the point Xq. To obtain the geodesic sphere we consider the set of solutions of the 
variational problem starting at point xq. We constructed the 4-dimensional projection of the 
geodesic sphere for the particle moving in the constant magnetic field. 



2 The equations of motion of a charged particle 

The optimization problem for the length functional on the set of the horizontal curves con- 
necting two given points is one of the main problems of the theory of optimal control [151 IS] ■ 
The velocity vector x'(t) of horizontal curves belong to the distribution A, hence we consider 
the metric tensor of the distribution A. If the vector fields ej(x), i = 0, . . .,3 form the basis 
of the distribution A(x), then its metric tensor is g%j{x) = (ei(x),ej(x)) x , i,j = 0, ...,3. 
Consider the length functional 

J(x(-),u{-))= [ L(x(t),u{t))dt, (12) 



to 



where 



/ 3 \V2 

L(x,u) = m | gij{x)u l u^ J (13) 

\i,j=0 J 

and m is the mass of the particle. Since the metric tensor has the signature (+, — , — , — ), 

the velocity vector ^j? belong to the cone of future V{x) = {v G A(x) \ (v,v) x > 0, v > 0}. 
This set depends of the point x. Hence let us consider another problem. 

All spaces A(x) are linearly isomorphic. Therefore we can consider the cone V(x) = {uG 

3 

M 4 | Yl 9ij( x ) uiv > J > 0) m° > 0}. For some neighbourhood of the point xq the inersection 

i,j=0 _ 

of all sets V(x) is nonempty and has a nonempty interior Q. In this neighbourhood of the 
point xq consider another problem for which u(t) G Q. If (x(t),u(t)) is an optimal process 



Note that nonholonomic variational problem is different from the nonholonomic mechanical problem. 
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for the problem (12), for which u(i) G V(x(t)), then this pair is locally an optimal process in 



the case u(t) G Q. From the horizontality condition (11) follows that for the velocity vector 
of the particle 

dx k 



dt 
dx A 
~dT 



u 



k = 0, 



3 



(14) 



^Aj{x)u>. 

3=0 



Let us consider the maximization problem for the length functional on the set of horizontal 



curves connecting two given points x l (to) = x l and x l (T) = x l T , 
- Pontrjagin function for the functional ((121, (14) has the form 



0. 



,4. The Hamilton 



H(x,u,p,a ) 



-a L(x,u) + 22 Pi 

j=0 



p*J2 a j< 

3=0 



X)U J 



(15) 



Due to Pontrjagin's maximum principle [TTl H~8] there is ao > and vector-function p k {t), 
k = 0, . . ., 4, to < t < T, satisfying the conjugate system of linear differential equations 



dpk 
dt 



dH 



dx k 



a 



(x(t),u(t),p(t),a Q ) 



dL 

dx h 



+ Pa 



3=0 



dx 



k = 0, 



(16) 



where (x(t),u(t)) is the optimal process for the problem (12). Since the metric tensor and 
4-potential of the electromagnetic field do not depend of the coordinate x 4 , the appropriate 
momentum component is the integral of motion (is conserved): p±(t) = const. From the 
equation ^ follows that p A is the length of the particle. For almost any t G [to, T) the 
function u \— > H(x(t),u,p(t),a ) reaches its upper bound on the set Q at u = u{t). If this 



maximum is reached at an inner point, then 



Pk - paM 



dH 

du k 

dL 



du 



k ■ 



0, i.e. 



k = 0, 



3. 



(17) 



Since the length functional does not depend on the parametrization of the path x(-), we 
can assume that the optimal control belong to the pseudosphere L(x,u) = 1. The normal 



covector for this pseudosphere has the coordinates 
impulse vector pk on the subset t "" ■• '' 



dL 



du h ' 

[u", . . ., u") has the coordinates pk 
function reaches its maximum when the projection of the vector pj, and the normal for the 



k — 0, . . ., 3. The projection of the 
p^Ak- The Hamilton 



unit pseudosphere are collinear. The equation (17) is exactly the condition for these two 



vector to be collinear. To exclude the vector p k from the equations (16) and (17) we use the 
identity 



A k (x(t)) - A k (x(t )) 



k = 0, 



3. 



Then 



a 



dL 

du k 



/ 7n^ dt ') +Pi S F i kV? dt ' + p*Mx(to)) - Pk(t ) = o. 

Jt Q OX J J tQ 



(19) 
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Parameters ao and Pk, k = 0, . . ., 4, can be multiplied by any positive constant. With ao = 1 
we obtain the equations of horizontal geodesies on the distribution A. If we additionnaly 
assume continuity of functions u(t), then 



d dL dL ^ . 

3=0 



dt du k dx k 



0,...,3. 



(20) 



These equations are identical with the equations of motion of a particle with the charge p± 
of the general relativity theory. 

For a = we obtain the equations of abnormal geodesies. Then at least one of the 
parameters Pk, k — 0, . . .,4, should be non-zero. Since in this case (17) p k = PiA k , k = 



0, . . ., 3, then p^ ^ 0. The solutions of the variational problem for ao = has the form 

3 



3=0 



0, k = 3. 



(21) 



The solutions of the variational problem (12) - (14) which do not satisfy Euler - Lagrange 



equations are called abnormal geodesies [T9l 120] . 

Montgomery constructed a distribution in 3-dimensional space for which the commutator 
of the basis vector fields is 2x 2 -^-% [211 E2J. In his example the abnormal geodesies are 

streight lines parallel with the x 2 axis. In our model the equations of abnormal geodesies 
on the distribution A has the form (21). If detF = 0, then this equation can posess a 
non-trivial solution. 



3 The geodesic sphere for the particle in the magnetic 
field 

The geodesic sphere has a compicated structure [21 [231 121] • In the following part of the 
paper we assume that M is the linear space M. 5 and the metric tensor of the distribution A 
is diagonal, g = Diag(l, —1, —1, —1). The equations of motion of a charged particle in the 



electromagnetic field (20) in the space with the diagonal metric tensor has the form 

m ^-^ F " M ' =0 ' fc = 0,...,3, (22) 

3=0 

where m is the mass of the particle, q is the charge of the particle, F is the tensor of the 

dx k 



dt 



electromagnetic field, u are components of the 4- velocity vector of the particle, 

u k , k = 0, ...,3. The fiveth component of the velocity vector is determined by the 
horizontality condition. 

Let us consider the motion of the particle in the magnetic field. Assume that the tension 
of the magnetic field H x ^ 0, and the tension of the electric field is zero. Then F23 = 
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Figure 1: Projection of the geodesic sphere on the coordinates (x 2 ,x 3 ,x 4 ). 



7^ 0, and other components of the tensor of the electromagnetic field F are zero. In 



this coordinates the equations of motion and the horizontality condition (11) has the form 

du 2 



m 



dt 
du 3 



qH x u A = 



m— h qH x u 2 = 

dt 



(23) 



dx 4 
~dt 



+ A ^ 



k=0 



Other components of the 4- velocity vector are constant. Assume that H x = const. Then 



from the equations (23) follows that in the magnetic field the charged particle moves along 
the helicoid line. 

Consider the projection of the geodesic sphere on the coordinates (x 2 ,x 3 ,x 4 ). In the 
plane (x 2 ,x 3 ) the particle moves along a circle. Choose the following components of the 
4-potential: A — 0, A\ = 0, A 2 = (px 3 , A 3 = 0. Then F 23 = — ip, H x = if. Designate 
P = rn'-P- Then the equations of motion of a particle in the subspace (x 2 , x 3 , x 4 ) has the form 



< 



' du 2 
dt 
du 3 
~dT 
dx 4 

, dt 



- pu 3 = 
+ pu 2 = 
+ (px 3 u 2 = 0. 



(24) 



The norm of the velocity vector along the geodesic is constant. Assume that (u 2 ) 2 +(u 3 ) 2 = 1. 
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Then 



x 2 (t) = — cos(a + pt) + b 2 



P 



x\t) 
At) 



sm(a + pt) + b 3 



P 

4 p 2 



(25) 



'-2pt + sin(2a + 2pt)) + -b 3 cos(a + pt) + 6 4 . 

p 



In the plane (x 2 ,x 3 ) the charged particle is moving along the circle with the center (62,^3) 
and the radius t— r. Let 6 2 = (cosa)/p and 63 = — (sina)/p. Choose the constant 64 so that 



\P\ 



x 4 (0) = 0. Then 



x 2 (t) 
x 3 (t) 
x\t) 



P 
1 



P 



(cos a — cos(a + pt)) 
sin(a + pt) — sin a) 



(26) 



^(— 2pt + sin(2a + Ipt) — sin 2a) + sina(cosa — cos(a + pt)). 



Ap 



The length of the geodesic 7 : [0, t] 



p- 

is equal to t. 



The equations (26) define mapping (t,a,p) t— > 7 QiP (t). In the plane (x 2 ,x 3 ) the point 
makes one turn around the center at \pt\ = 2n. More long geodesies (\pt\ > 2n) are not 
the shortest paths. Here one should consider the shortest paths, since the restriction of the 
metric tensor on the tangent bundle for the coordinate plane (x 2 , x 3 , x 4 ) is positively defined. 
The geodesic sphere of the radius s can be defined by the following formula: 



B(s) = {l a , p (t) I a e [-it, it], t = s, \pt\ < 2tt}. 



(27) 



Fig. 1 shows the projection of the geodesic sphere of the radius s = 1 for the coordinates 
(x 2 ,x 3 ,x A ) at ip — 1, with p e [— 2tt,0] on the left and p G [0, 27r] on the right. If \pt\ < 2tt, 
the geodesies are the shortest paths. The nonholonomic geodesic sphere of any radius has 
a pair of points such that if we continue the geodesic after this point we enter the geodesic 
sphere with some lesser radius. If \pt\ > 2rr, then the geodesic of length t ceases to be 
optimal path. All possible endpoints 7«, p (t) of the geodesies of length t reach the x A axis 
at pt = 2irk, k e Z, this point being the intersection point for the considered surface. The 
density of these points increases when we approach the origin. This phenomena is called the 
nonholonomic wave front (fig. 2). Fig. 2 shows the projection of the surface defined by the 
geodesies of length s = 1 for the coordinates (x 2 , x 3 , x A ) at ip = 1 with p £ [— 87r, — n] on the 
left and p £ [tt, 811} on the right. These geodesies cease to be optimal paths for \p\ > 2n. For 
large enough p this surface belongs to some cone with the apex at the origin. This cone has 
the axis of symmetry x A . Indeed for p large enough from the equations (26) follows that 



At) 



2p 



(p 



00 



(28) 



In the plane (x 2 , x 3 ) the particle moves along the circle of the radius r^-. The distance from 



the starting point to the ending point is changing from to 



\P\ 



For p large enough the 
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Figure 2: The non-holonomic wave front defined by non-optimal geodesies. 

geodesies of the length s belong to the interior of the cone with the angle of inclination ^-s 
and the apex at the origin. 



4 Distance along the x A coordinate 

In the subspace (x 2 ,x 3 ) the particle moves along the circle. Let us find the values of t for 



which the end of the geodesic is over the starting point in the subspace 



x , x 



, x ). The end 



of the geodesic is over the starting point for pt 
the equation (26) at the point p = + 



■ 2nk, k G Z. Use the Taylor expansion for 
0: 



x\t) 



i 1 



t 



-ip 



Aixk 



+ tP ^ 9 + 0{6) ( ^ 0) ' 



(29) 



Since t = s, the distance for which the particle moved along the x 4 coordinate is of the order 

^iTfk' 

For distribution with the positively defined metric tensor there is the theorem about 



the ball shape [22J [26]. In the coordinate neighbourhood of any point one can choose the 
vector fields e% which form the basis of the distribution A. The field of planes that span 
commutators of all possible vector fields of the distribution A is designated Az = [A, A]. 
One also consider Ak+i = [Ak,A]. This sequence is stabilizing: there is the minimum m 
such that A m = A m+ i. The number m which can depend of the point of the manifold is 
called the degree of nonholonomity of the distribution A. Designate = dim Ak and define 
the function ip(i) = j, if rij-i < i < rij. The function ip is defined on the set of numbers 
{1, . . ., dim^4 m } with values at {l,...,m}. For completely nonholonomic distribution on 
the Riemannian manifold there are the coordinates x % such that the ball of accessibility 
is bounded from above and from below by the set \x % \ < e^. For the two-dimensional 
distribution on the 3-dimensional manifold Uq = 0, n\ = 2, n 2 = 3. Therefore the function 
ip{l) = 1, (p(2) = 1, <p(3) = 2. 
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5 Lagrange formulation for the equations of motion 

5.1 General case 

Let us consider the classical problem of the calculus of variations: how to find an absolutely 
continuous vector-function which maximizes the functional [2J [T7] 

J(x(-),u(-))= [ L(t,x(t),u(t))dt, (30) 



where u(t) = tt. Assume that x(0) = Xq, x(T) = x±, and allowed paths satisfy the 

conditions (f l (t,x, ^) = 0, i — 1, . . .k. Then there are k measurable and limited functions 
Aj(t) called Lagrange multipliers and a constant do > 0, at least one non-zero and such that 
the function 

k 

L x (t, x(t),x(t)) = a Q L(t, x(t), x(t)) + Ai(t)^(t, x(t),x(t)) (31) 

i=i 

almost anywhere on [0, T] satisfies the Euler - Lagrange equations in the integral form 
dL x (t, x(t),x(t)) f T dL x (r, x(t), x(t)) 



dx\t) 



dx 1 



dr + Q, i = l,.. .n, (32) 



where c, are constants [TTJ p. 279]. The parameters a , Aj can be multiplied by any positive 
constant. With a = 1 we obtain regular geodesies. With a = we obtain abnormal 
geodesies. 

Now assume (for the general case) that the distribution B is defined by the family of 
differential forms lo\ j = m+1, . . .n. The horizontality conditions have the form 

ct^(u) = 0, j = m+1, . . .n, x = u. (33) 

Then the Lagrange function for the considered problem is 

n 

L x (x,u) = a (u,u) 1/2 + Y W(w), (34) 

j=m+l 

where (■, •) is some non-degenerate^] bilinear form on TM n . The equations of the horizontal 
geodesies are 

Du n ■ n 

a °^'')+ Yl X ^ J+ W(«:0=O. (35) 

j=m+l j=m+l 

where is the covariant derivative of the velosity vector of the particle along the path. 
The velocity vector can be decomposed by the basis vector fields of the distribution: 

u(t) = J2v k (t)U<t))- (36) 

k=l 



3 This assumption for our model is not necessary, as shown below. It is enough for the restriction of the 
metric tensor on the distribution to be non-degenerate. 
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Then the covariant derivative is 

^£^+5>^, (37) 

k=l i,j=l 

where V is the symmetric connection defined by the bilinear form. Use the Maurer - Cartan 
formula 

dw(^v) = &(t})-Vu(.Z)-u([S,v])- (38) 

Since u^(^i) = 0, % = 1, . . .m, j = m+1, . . .n, then for the velocity vector and the basis 
of the distribution cLu^u,^) = — uji([u, &]). This commutator can be rewritten using the 
structural constants of the distribution B: 

m m n 

[u, &] = Yl v ^i] = X> £ c k&> ( 39 ) 

k=l k=l s=l 

n 

where = ^2 c ti£s- This basis and the differential forms can be selected such that 

3=1 

o; J (£j) = 5{, i,j = m+1, . . .n. Then for the distribution B 

a °<^>&>+ E A iE^ = °. Z = l,...m. (40) 

j=m+l k=l 

Since ^(^i) = Sj, i,j = m+1,. . .n, then for the velocity vector and all basis vector fields 
of the tangent space du^(u,^i) = — a;- 7 ([it, &]). Hence we can write the projection of the 



equation (35) at the basis vectors which do not belong to the distribution: 
,Du 

j=m+l k=l 



a ( — ,£i) + \i+ A jE^ = ' l = m+l,...n. (41) 



5.2 Lagrange formulation for the distribution A in our model 

3 

The distribution A is defined by the differential form u x — Aj(x)c?x J + c?x 4 . The Euler 

i=0 

Lagrange equations for the length functional (12) with the condition u^^'it)) = are 



a (^,-) + Acu + A^(7,-) = 0, (42) 
where the covariant derivative 

% = ±%* + ±<"^«. («) 

fc=0 ij'=0 

V is the symmetric connection defined by the metric tensor. Please note that in the Lagrange 
method we have to use the metric tensor defined on the whole tangent space TM whereas in 
the maximum principle we used the restriction of the metric tensor on the distribution only. 
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This is no problem because in the result we use the components of the metric tensor of the 
distribution only. We can also use the connection (Christoffel symbols) with all lower indexes 
and in this case the bilinear form does not have to be non-degenerate. The velocity vector 

3 

7(*) = E v k (t)e k (-f(t)). The Euler - La grange method also has the constant ao > [T7] 

k=0 

p. 279]. The parameters a , A can be multiplied by any positive constant. With a = 1 we 
obtain the equations of regular geodesies. With ao = we obtain the equations of abnormal 
geodesies. Projecting the equations of motion on the basis of the distribution A we obtain 

D' 3 

a (-^, + 4 vj = °> * = 0, . . 3. (44) 

j=0 

4 

The structural constants are defined by = c \j£ii where £j = ej, j = 0, . . ., 3, £ 4 = <9 4 . 

/ o 

For the distribution A [e*, ej] = —F^d^, [e 3 -, d&] = 0. Therefore cfj = Fji, and other structural 
constants are zero. Since (x> x (<9 4 ) = 1, then 

a (^,d 4 ) + X = 0. (45) 
at 

We can assume that (^c^) = 0. Then A = const, and A can be interpreted as the charge 
of the particle (or charge to mass ratio depending on the Lagrangian). Note that applying 
the maximum principle in section [2] we used the restriction of the metric tensor on the 
distribution only. Since this restriction is non-degenerate, we can rise appropriate indexes 
of the Christoffel sybols and get ordinary differential equations. The result does not depend 
on the extension of the metric tensor on the whole tangent space. 

All figures published in this paper were produced using our own 3D graphics program 
©V.R. Krym. 
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